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1 INTRODUCTION 



ABSTRACT 

In 2004, McLaughlin et al. discovered a phenomenon in the radio emission of PSR J0737— 3039B 
(B) that resembles drifting sub-pulses. The repeat rate of the sub-pulses is equal to 
the spin frequency of PSR J0737— 3039A (A); this led to the suggestion that they 
are caused by incidence upon B's magnetosphere of electromagnetic radiation from A. 
Here we describe a geometrical model which predicts the delay of B's sub-pulses rela- 
tive to A's radio pulses. We show that measuring these delays is equivalent to tracking 
A's rotation from the point of view of an hypothetical observer located near B. This 
has three main astrophysical applications: (a) to determine the sense of rotation of 
A relative to its orbital plane; (b) to estimate where in B's magnetosphere the radio 
sub-pulses are modulated and (c) to provide an independent estimate of the mass ratio 
of A and B. The latter might improve existing tests of gravitational theories using this 
system. 
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The discovery of the double pulsar system PSR J0737-3039 (IBurgav et al. 20031 : iLvne et al. 20041 ) has led to important 
advances in the study of radio pulsars and their associated phenomenology. The independent determin ation of the orbit o f 
both neutron stars, and the recent measuremen t of five post-Keple rian parameters from the timing of A ( Kramer et al. 20061 ) 
and a sixth one from eclipse observations of A I Breton et al. 200^ has made this the most over-constrained system known, 
allowing five tests of general relativity. 

PSR J0737— 3039B (henceforth B) has many unique f eatures. It emits pulsed radio waves during most of the orbit, but 
it becomes much brighter at two distinct orbital ph ases (|Lvne et al. 2004 ). Furthermore, its pulse profile changes strongly 
with orbital phase and with time ( Burgav et al. 20051 ). Though instructive in other ways, this behaviour has made it difficult 
to measure times of arrival accurately and obtain a precise estimate of the semi-major axis of the orbit of B, limiting the 
precision of our knowledge of the pulsar mass ratio R = mA/mB and the precision of some of the tests of general relativity 
in this system (see, e.g., Kramer & Stairs 2008 for a recent review). 

This paper is motivated by anoth er unique feature of B: fo r certain orbital phases its radio emission is clearly modulated 
by electromagnetic emission from A (jMcLaughlin et al. 20041 ): the phenomenon superficially resembles drifting sub-pulses 
(e.g. Lorimer & Kramer 2005). In this work, we discuss the timing of this phenomenon. In § [2] we present a theoretical 
calculation of the delay in the arrival at the Earth of the sub-pulses of B relative to the radio pulses of A responsible for the 
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Figure 1. Sequence of events in A, looking down the pulsar's spin axis. We assume that the plane of the drawing is the orbital plane. 
Left: at to the radio beam is pointing at the Earth (outside of the orbital plane), producing a maximum of radio emission. At this time, 
A's EM beam is not necessarily pointing at B. Center, top: at ij, A is emitting the EM signal that will produce the response. Locally, 
the wavefront (light line) is a plane perpendicular to the direction of propagation of the EM signal. Between to and t'^, A rotated by 
9 = Aj, + 7r/2 — 0e if the pulsar is rotating in the same sense of the orbit (clockwise here), or 6 = 2tt — (A;, + it/2 — (f)^) if the pulsar is 
rotating in the opposite sense of the orbit. Right, bottom at time t'-^, the EM wavefront arrives at B's position. Soon after that, at time 
t\ , the wavefront arrives at point E, where the Earth-bound radio emission of B originates. This is located at a fixed distance e from B 
on that pulsar's line of sight to Earth, it is at a small distance from the orbital plane ecosi (see text). The arrival of the wavefront at 
point E at time ti produces the response in B's emission, which then follows to Earth. Aa is the longitude of A as seen from E at ti ; so 
•^6 = -^a + TT is the longitude at which A's EM beam transmits towards B. 



drift. We name this simply the "response delay". In §[31 we discuss how we can compare the predicted and measured response 
delays. In § O we highlight the astrophysical knowledge that can be gained from the timing of the response delays. 



2 CALCULATING RESPONSE DELAYS 

2.1 Time between emission of A's radio pulse and B's response 

The time between the emission of A's radio pulse and B's response can be divided in three parts, which we discuss below. 



(i) The sub-pulses of B are spaced by 22.7 ms, which is the same as the spin period of A (jMcLaughlin et al. 2004 ). Therefore, 
we can think of the maximum intensity of the sub-pulses of B (the "response") as being caused by the impact upon B's radio 
emitting region of a co-rotating beam of electromagnetic radiation from A. We henceforth refer to this as the EM beam. Since 
the EM beam is not necessarily aligned in rotational phase with A's radio beam, as shown in Fig. [1] we introduce a phase 
offset <j)e. between the two beams. We define A's zero longitude to be the point where the radio emission towards the Earth 
reaches maximum intensity. As we show later, one of the applications of our model would be the direct determination of cj>e. 

One of t he simplifying assumpti ons in Fig. [T]is that the spin axis of A is perpendicular to the orbital plane. This has been 



suggested ([Manchester et al. 20051 ) as the most plausible explanation for the lack of change in the pulse profile of A; Ferdman 
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(2008) set an upper limit of about 14° for the misalignment between the orbital angular momentum and the spin of A. We 
henceforth assume that the spin axis of A is indeed perpendicular to the orbital plane. It is not know whether A rotates in the 



same sense as the orbit or in the opposite sense, but kinematic considerations ( Bailes 19881 ) make the former more likely. We 
discuss this issue further in § |3] and § |4l 

In Fig. [T] at time to we see A's radio beam pointing at the Earth, causing the radio emission to reach a maximum. At time 
to the EM beam is pointing at a direction At, where it emits the signal that will later cause B's respons^. Between the two 
events A has rotated by an angle 6 = ±(Aa — n/2 — 4>e) where Aa = At — tt is the longitude of A as seen from the location and 
time where the EM signal modulated the radio pulse of B (point E at time t\, see Fig. [l]). The positive sign corresponds to a 
clockwise rotation of A as shown in Fig. [T] i.e., in the same sense of the orbit. The time between the two events is given by 

t'o-to = PA^, (1) 

where Pa is the spin period of A. 

(ii) After t'g the EM pulse of A travels towards B; at time t'l it gets to B's exact position. The distance between these two 
events is r^s(Aa), this is calculated in ij 12.31 In the assumption above (that signal travels at the speed of light) we have 

t[-t',= '-^^. (2) 

(iii) What happens near B is less certain. The simplest assumption is that the EM signal is able to travel through B's 
magnetosphere at the speed of light (this is possible e.g. if the EM signal consists of high-energy photons) and modulates B's 
Earth-bound radio emission at the point where it is being generated; we designate this as "E". 

We must keep in mind that this simplified assumption is not the only possible case: a) the modulation could occur after B's 
radio signal is produced, i.e., somewhere between point E and the Earth and b) more generally, the chain of events from the 
emission of the EM signal at A (tg) to the modulation of the radio signal of B (ti) might not lie along a single null geodesic. 
For the time being, we will quantify the simplest case only (see Fig. [1} . 

Since E is where the Earth-bound radio emission is being generated, we assume it is the point of the radio-emitting region 
of B's magnetosphere that is in the line of sight from B to Earth. This only happens when B's radio pulse is "on". For most 
of B's rotational cycle, when its pulse is "off", no part of its radio-emitting region is pointing at Earth, and no response is 
detected. 

In our simplified model, we assume that the distance of E to B is a constant e as a function of B's spin phase. Therefore, in 
the reference frame of the orbital plane (with axes along the line of nodes, perpendicular to the line of nodes and perpendicular 
to the orbital plane), E is located at fixed coordinates (0, esini, ecosi) relative to B, where i is the orbital inclination of the 
system (see Fig. [TJ . 

At t'l, the EM wavefront from A reaches B. A small amount of time later, at ti, the same wavefront reaches E and produces 
the response. From that Figure, we can see that 

ti — t'l = -sinAaSini. (3) 

c 

Between t'l and ti the orbital motion of B will make E move by a small amount relative to the center of the binary, this will affec t 



the arrival time at E. We note, however, that for slow pulsars e/c is expected to be of the order of 1 ms (jBlaskiewicz et al. 1991h . 
Assuming this is true, a detailed calculation of the 0(ec~^) Doppler correction due to B's (and E's) motion shows that this 
term is smaller than 1 fis. Since the r.m.s. timing precision of A is about 10 jUS, and the timing precision for B's response is 
likely to be worse, we will ignore all terms with magnitude < 10 fis. 



For the total time difference we obtain: 



ti-to = {ti -t'i) + (t'l -t',) + {t'a~ to) = - sin A„ sini + li^^lM + + 0{ec-^), (4) 

which should be valid in any inertial reference frame using the orbital parameters as measured in that frame. 



2.2 Response delay 

We now calculate the delay between the reception at the SSB of A's radio pulses and B's responses, A(Aa). To do this, we will 
start by adding to £i — to the difference in the ranged at which these events occur, zi — zo divided by c; this is the so-called 
"R0mer delay". 

^ Lower-case subscripts refer to A computed for the emission-reception-Earth triangle, while uppercase subscripts refer to lambda 
computed for a line through the centre of mass of the binary system, O. Here we have neglected aberration effects and the small angle 
between AB and AE, as the corresponding time delays are below the accuracy relevant for this paper. 

We use the term "range" to indicate the distance of a given event from the center of mass of the binary projected along the direction 
to the SSB as seen at the SSB. 
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Figure 2. Emission (at tg, on A) and reception at the position of B (at t'^) of the EM signal that causes B's responses. Upper diagram: 
events depicted in the orbital plane. O represents the centre of mass of the system, orbital motion is represented as clockwise. Lower 
diagram: events depicted on a plane that is perpendicular to the orbital plane and the plane of the sky. This contains the line of sight 
to the Earth and allows a representation of the response delay s. The orbital inclination i of the PSR J0737— 3039 binary system is much 
closer to 90° than depicted in this figure l lKramer et al. 2006h : we represent it as being substantially lower for clarity. The difference in 
the locations of to,tQ and is too small to be discerned at this scale. In this figure, the line of periastron just happens to coincide 

with the line of nodes, with u)j\ = 180°, = 0°. 



As in the case of t\ — to (eq.|4ll, zi — zo can also be described as the sum of three terms, {zi — z'l) + {z'l — zq) + {z'q ~ zq): 

(i) The difference in range between the events at to and tg is given by: 

z'o~Zo = {t'o~to)vi^A=PA—Vl,A, (5) 

where Vi^a is the velocity of A relative to the centre of mass, projected along the line of sight to Earth ( H2.3|l . 

(ii) In the lower diagram of Fig. [51 we see that the difference in range between the events happening at t'-^ and t'o is given 
by: 

z'l — z[) = —rAB{\a)sttl\aSVai. (6) 

(iii) Ignoring the motion of E between t'^ and ti , we have 

zi-z[ = -e. (7) 
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The delay is therefore given by: 

Ar{Xa) = ti ~ to -\ = (1 -smAaSmi) +Pa7; h 0(ec ), (8) 

c c 2n 

where the subscript "r" indicates that we are taking into account the R0mer delay only. In this equation we use P4 — 

Pa{1 + vi^a/c) to represent the spin period of A with the classical Doppler correction due to its orbital velocity. This equation 

is valid at the SSB if we use the orbital parameters as measured there. 

In Appendix A we present a detailed calculation of the relativistic contribution to the total response delay A, known as 

the "Shapiro delay" (As). For the range of orbital phases where we observe B's responses As changes 4/^s, i.e., smaller than 

the r.m.s. timing precision of A. We can therefore assume A = Ar + As — Ar for the remainder of this paper. 



2.3 Distance between A at transmission and B at response 

Having calculated the response delay as a function of Xa and the separation between the two pulsars, rAs(Aa), we now 
calculate these two quantities as a function of Xa at t^. The reason for this is that Xa can be calculated precisely for any 
given time using the equations in Damour & Deruelle (1985,1986), furthermore, this quantity determines the instantaneous 
configuration of the system. 

In the equations that follow, we will only use the Newtonian terms to 0(c~^). A detailed calculation shows that the 
0(c~^) Newtonian terms are much smaller than the r.m.s. timing precision of A. This means that they can be ignored for or 
present purposes. 

At t'l), the separation of pulsar I (A or B) from the centre of mass is given (see, e.g., Roy 1988) by 

l-e" 

di = ai— -, (9) 

1 + e cos / 

where aj is the semi-major axis of its orbit, e is the orbital eccentricity and the angle / is the true anomaly of the system at 
^0 if ~ Xi — LOi , where uji is the longitude of periastron of pulsar I at t'o). T he compon ents of a pulsar's velocity along the 
radial (i.e., away from the centre of mass) and transverse directions are given ( Rov 198^ ) by 

di — w/esin/, (10) 

dif = «/(l + ecos/), (11) 

where 
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2-K ai 

is the transverse velocity of pulsar I at quadrature (/ = 7r/2) and P(, is the orbital period of the binary. The velocity of that 
pulsar relative to Earth is given by (e.g., Green 1985): 

vij = VI sin i (cos A/ + e cos uoi) . (13) 

Calculating this velocity for A we can immediately quantify the Doppler correction term for Pa in eq. |8l 

To calculate the instantaneous distance between the pulsars (cJas), we replace a/ by a = + as in eq. |9l During the time 
it takes A's signal to cross this distance (dyis/c), B is moving. Because the EM signal is moving in the radial direction, only the 
radial component of B's motion (dr) will affect the travel time between the two pulsars. During the cross time dr ~ dBdAB/c, 
i.e., the distance between the events at and t'l is 

TABiXA) = dAB +dr = dAB (l + — esiu/) + O (c-2) ; (14) 



the perpendicular motion of B between smd t'l is 

dr± = dBdf ^ dBfit'i - to) = «s(l + ecos /■)!ji£(M. (15) 

c 

Looking at Fig. Owe can see that 

Xa = \A+ , =\A + —il + ecosf). (16) 

rAB(-\A) c 

In the derivation above we made the approximation that Xa is the longitude of A as seen from B at t'l. The longitude of A as 
seen from point E at ti is larger by a very small amount, esin i cos Aa/rAB- 
Noting that dx = dr± cos Xa , 

rAs(AA) sinAa = TABiXA) siu Aa + dr_i_ cosAa. (17) 
With this result we can re- write eq. [8]as a function of Aa: 

A(Aa) = '^^B{XA)-e (^_gjjj^_^gi^ •) _ ^cosAAsini + P^;^ + 0(ec-'). (18) 

C C ZTT 

Using eqs. [T^ and we can re- write this as a function of a Keplerian term, K{Xa), which can be calculated from known 
orbital parameters and the system's geometry, and the unknown quantities e and 9: 

A(Aa) = Jf(AA) (1 - sin Aa sin i) + Pa 7^ -I- 0(ec-^) (19) 



K{Xa) = -(1 
c 



1 — sin Aa sin i 
1 + e cos / 



( ^ vb . ,\ Vb \ ■ • 

1 H e sm i cos Aa sm i 

V c J c 



+ O (c-^) . (20) 



The K{Xa) term is the time difference, measured at the SSB, between the events that occurred at t'l and to in the reference 
frame of the binary; it is by far the largest contribution to the response delay. 



3 COMPARING MEASUREMENTS OF RESPONSE DELAYS WITH PREDICTIONS 

So far, we have just made a theoretical calculation of the response delays, A(Aa). We now discuss what we can learn from 
actually measuring such delays. 



3.1 Measuring absolute delays 

If we measure A(Aa) for several Aa in a single orbit and subtract K[Xa), we can measure the delays that are not a priori 
predictable: 

A(Aa) - ^(Aa) = - (sin A„ sini - 1) ± PA f ^ - 1^ - iV (21) 
c \27r ZTT 4/ 

In order to measure A(Aa), we have to specify a particular radio pulse of A (in principle the closest to the time of emission of 
the EM signal that caused that particular response) and subtract its (barycentric) time from the response's barycentric time. 
However, because we have no prior knowledge of the terms in eq. 1211 there is some ambiguity in the choice of the pulse that 
is closest to that time of emission. 

An experimental method that does away with the issue of the choice of a pulse of A and even the measurement of A(Aa) 
entirely would proceed along the following lines: 
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(i) B's responses are p articularly noticeable in an intensity gray-scale plot of time versus spin phase of B as presented by 
I McLaughlin et al. 2004 ') . In such a plot, we can determine the precise spin phase corresponding to each response; the precise 
method to accomplish this and the attainable precision will be discussed elsewhere. 

(ii) In the next step, we convert this phase into a barycentric time. This can be achieved using the software package tempcH. 
The result corresponds to an event that occurred at ti in the reference frame of the binary. 

(iii) Subtracting K{Xa), we obtain a second barycentric time that nearly corresponds to A's emission of the EM pulse 
(to + e/csin Aa sini) in the reference frame of the binary). At this time we calculate Xa from the binary's ephemeris. 

(iv) Subtract the nearest time at which the rotational phase of A is zero. This time should be the equivalent of A(Aa) — 
K{Xa)- By fitting this to eq. [2T]we should be able to determine e, 0^ and the direction of the spin of A. 

With data from different orbits, we can also get a sense for the stability of and e. Thus far we have assumed that e 
(the height above B where the response is produced) is a constant as a function of B's rotational phase. It is possible, however, 
that e varies with B's spin phase. This would also introduce a secular variation of the e observed for any spin phase of B due 
to the pulsar's geodetic precession. If that is the case, we could map the height of E above B as a function of the precession 
phase. Confusion with K{Xa) can be avoided because, despite the fact that both terms vary with 1 — sinAAsini, only the 
latter varies with (1 + ecos/)"^. 



3.2 Delay variations at a constant Xa 

The response delays A(Aa) will vary in time for any given Xa because of the apsidal motion ui changes /, which causes a 
change in K{Xa)'- 

6t^-t,{^A) = At=t,(XA) - At=ii(AA) = -(1 - e") (1 - sinAAsini) [j(t2) ~ j{ti)], (22) 

c 

where ti and t2 are two epochs at which the longitude Aa occurs and 

j{t) = (l + l^esin/) (l + ecos/)-\ (23) 

Other terms of order larger than cancel out in the subtraction because of the constant Aa, and the same happens for the 
terms dependent on 6 and The variation of the Doppler correction to Pa are < 1 /is and can therefore be ignored. 

The St2-ti (Aa) are important because, being due solely to a variation of K{Xa), they are proportional to a. The maximum 
variation occurs between the time when Aa coincides with periastron (/ — 0, occurring at t — tp) and apastron (/ — n, 
occurring at t — ta, which for PSR J0737— 3039 is about 10.65 years later): 
2ae 

(5t,-tp(AA) = — (1 - sinAAsini) . (24) 

In a single 2003 observation, McLaughlin et al. (2004) observe B's responses from ~ 195° < Ab <~ 225°, i.e., occurring 
shortly before the magnetosphere of B eclipses A at As = 270°. This corresponds to ~ 15° < Aa <~ 45°. At these extremes 
we have St^-tj,{XA) ~ 0.3814(3)s and 0.1508(l)s. The uncertainty in this prediction is entirely due to the uncertainty of a 
(which results entirely from the uncertainty in the measurement of as); it represents 0.5 - 1.5% of a rotation of A; this means 
that we should be able to keep track of the response delays with very high confidence. 

The variations themselves are certainly measurable, since they are equivalent to many rotations of A. This means that 
we can always make an independent measurement of a. Depending on the precision and number of measurements of A(Aa), 
the value of a derived in this fashion might be more precise than the present value derived from timing. In that case, we can 
improve our knowledge of the mass ratio: 

R=^-l, (25) 

CXA 

where xa is the projected semi-major axis of A's orbit, in light seconds; this quantity is directly and very precisely measured 
from the timing of A. This equation (and also eqs. I22l andl24 p show that, to determine ii from a and xa we need to know sini. 
We could in principle determine sini independently from comparing measurements of 5t2-ti{XA) made at different orbita l 



longitudes . However, most theories of gravitation predict the Shapiro "s" term to be the same as sin i (|Damour &i Taylor 1992 



Will 19931 ): therefore such an independent determination is not likely to be a useful test of gravitation. For this reason, we 



can use the extremely precise estimate of sini from s in the estimate of R. 

Finally, we remark that there is a small difi'erence between the intrinsic and observed semi-major axes of the pulsar's 
orbits that is caused by aberration effects. Kramer et al. (2006) estimate that for A these are of the order of 10~^xa, i.e.. 



^ http:/ /www. atnf.csiro.au/research/pulsar/tempo/ 

■* The latter term is likely to cancel even if the main assumption in il2.1l is not correct, i.e., it only requires that any terms dependent on 
e vary with Aa in a repeatable way. 



8 Freire et al. 



about 1/is; this is similar to the level of precision for the measurement of xa- It is highly unlikely that we will measure a to 
a similar precision, so this should not be a problem for the determination of R in eq. 1251 For B, Kramer et al. estimate that 
the effect of aberration should be of the order of 1Q~'^xb- It might be possible that we measure xb = asini/c — xa to better 
than the necessary level of precision. However, to be able to measure the aberration directly, we must be able to measure 
Xb from timing with this level of accuracy, so that the differ ence between the t wo values becomes evident. The intrinsic 
uncertainty of the xb obtained from timing is about 10~^a;s ( Kramer et al. 2006h . one order of magnitude larger than the 
effect of aberration. Therefore, unless the direct timing of B improves by more than one order of magnitude, the effects of 
aberration won't be separately measurable. 



4 IMPLICATIONS AND PROSPECTS 

In this paper we have calculated the delay between the radio pulses of A and the modulated radio pulses of B, A(Aa), assuming 
a simple scenario for that modulation, i.e., that it happens at the point, E, where the radio emission of B is being produced 
and that the modulating signal is traveling from A to E at the speed of light. Things could be more complicated, particularly 
considering that this trajectory must intersect the magnetosphere of B. 

We will present details of the measurements of A(A,4) and their timing analysis elsewhere. If the response delays obey 
the equations presented above in a consistent manner, that will validate our model. In that case, the measurement of these 
delays relative to the radio pulses of A can provide us with new astrophysical information. 

First, a determination of the sense of A's rotation relative to its orbit, something never achieved for any other pulsar, 
would introduce fundamental constraints on binary evolution scenarios for this pulsar, as well as improved constraints on B's 
supernova kick. Findi ng that the o rbital angular momentum is anti-aligned with A's angular momentum would require a very 



large supernova kick (|Bailes 1988f ). Previous studies (e.g., Willems et al. 2006, Stairs et al. 2006) predict instead that the 



kick that produced B was rather small, therefore the angular momenta should be aligned. A confirmation of this alignment, 
combined with the small angle between the momenta, would introduce stringent constraints o n the magnitude of the k ick. 
Furthermore, we would know the sign of the expected relativistic spin- orbit contribution to a; ( Damour fc Schafer 198^ ^ It 



is probable that this contribution will be measured in the near future jKramer fc Wex 20091 ): that would allow an estimate 



of A's moment of inertia. The moment of inertia of A, together wi th the well-determined m ass of the pulsar, will introduce 
fundamental constraints on the equation of state for dense matter ( Lattimer &: Schutz 2005h . 



Second, by measuring e, or by introducing upper limits to it, we might be able to locate the region where the EM signal 
from A is modulating the radio signal from B and relate it to the region where we expect the radio emission of B is being 
generated. Thus far, the best location of pulsar radio emission comes from the interpret ation of multi-frequency polarimetric 



pulse profiles in light of a relativistic version of the the familiar rotating vector model (|Blaskiewicz et al. 19911 ): the results 
indicate an emission height of a few hundred km. If the signal is being modulated as its being generated, then e should be of 
the order of a few ms and therefore a measurable quantity. If e is significantly larger, then the modulation is happening after 
B's radio signal is generated. 

Third, we will make an independent measurement of the orbital separation of the two pulsars, a. If the timing of B's 
responses is precise enough, this might give us a more precise measurement of the mass ratio of A and B, which would increase 
the precision of some of the previous tests of general relativity carried out in this binary system. This is a key input parameter 
for the general tests on conservative gravity theories outlined in Kramer & Wex (2009). 

In the ideal case that we can track the response times well, one might think of their reception at the Earth as being 
equivalent to having a radio telescope at an altitude e above B tracking A's rotation and then relaying the results live to Earth. 
This analogy highlights how fortunate we are to have a phenomenon like B's responses. 
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APPENDIX A 



Including post-Newtonian corrections, the time taken for a photon to propagate from the point of emission xo to x in the 
gravitational field of a TV-body system is given by 



t-to = 



X - Xo 



2gA , 

1=1 



rox + ri + R 



roi +ri - R 



(26) 



where rrn is the mass of the i-th body located at Xi, roi = |xi — xo|, rt = |x — Xi|, and = |x — xoj (jBrumberg 199lf l. The 
first term, |x — xo|/c, is the R0mer term used in the main text. The second term, the relativistic correction, is the sum of the 
Shapiro delays caused by the individual bodies in the A^-body system. 

We now apply equation (|26|) to the signal propagation discussed in this paper (see figures [T] and [Sjl . Neglecting terms 
smaller than a few fis one finds for the signal propagating from A to Earth at distance D 



,(A®) 



2GmA 



In 



r_D 



2Gmi 



In 



2D 



dAB{l + cos A 



(27) 



Sr. (^ (Iab) is the emission height of the radio signal in the magnetosphere of A and Aa,e is the angle between the direction 
to Earth and the direction to pulsar A as seen from B. For the signal propagating from A to the point E to sufficient accuracy 



.(AE) 



2GmA 



In 



dAi 



2Gms , 
H ; — m 



2dAL 



e(l + cos \a.( 



(28) 



5em (^ d,4s) is the emission height of the EM signal that triggers the sub-pulse emission at E. Finally, for the sub-pulse 
signal emitted at E the Shapiro delay is given by 



2GmA 



In 



2D 



2GmB , 
H T, — In 



dAs(l - cos Aa,®) 
The contribution of the Shapiro delays to A is therefore given by 



2Gm, 



In 



1 -I- sin i sin Aa 



, 2GmA , 
H 5 — In 



(5em 



AGniB , 
H 5 — In 



dAi 



(29) 



(30) 



where we have used cosAa.e ~ — sinisinAA. The second term of equation (|30p is constant and, to sufficient accuracy, can 
be absorbed in <j!>e. The third term of equation (|30p changes along the eccentric orbit with an amplitude of just 2 fis, and 
therefore can be absorbed into (p^ as we ll. In principle th e first term of equation (|30p can change quite significantly along the 
orbit, as sini ~ 1 in the double pulsar ({Kramer et al. 2006. 1. However, as mentioned above, the sub-pulses are only observed 
in the range ~ 15° < Aa < ~ 45°. Across this interval the first term of equation (|30|l changes by less than 4 ^s. Consequently, 
Shapiro delays can be ignored in the calculations of this paper. 



